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ABSTRACT 

We present f/ST/FOS 1 spectra of the two bright images (A and B) of the gravitationally lensed QSO 
0957 + 561 in the wavelength range 2200-3300 A. We find that the absorption system (z abs = 1.3911) near 
z em is a weak, damped Lya system with strong Lya absorption lines seen in both images. However, the 
H i column densities are different, with the line of sight to image A intersecting a larger column density. 
The continuum shapes of the two spectra differ in the sense that the flux level of image A increases more 
slowly toward shorter wavelengths than that of image B. We explain this as the result of differential 
reddening by dust grains in the damped Lya absorber. A direct outcome of this explanation is a determi- 
nation of the dust-to-gas ratio, /c, in the damped Lya system. We derive k = 0.55 + 0.18 for a simple 1/A 
extinction law and k = 0.31 +0.10 for the Galactic extinction curve. For gravitationally lensed systems 
with damped Lya absorbers, our method is a powerful tool for determining the values and dispersion of 
k , and the shapes of extinction curves, especially in the FUV and EUV regions. We compare our results 
with previous work. 

Subject headings: galaxies: ISM — gravitational lensing — quasars: absorption lines — 
quasars: individual (QSO 0957 + 561) 


1. INTRODUCTION 

0957 + 561 A, B, discovered by Walsh, Carswell, & 
Weymann (1979; hereafter WCW) are the first examples of 
multiple images of a gravitationally lensed quasar. This 
system has an emission redshift of 1.4136 (Weymann et al. 
1979), and the separation between the two images is 6" 1 7 
(Roberts et al. 1985). WCW noticed that below 5300 A the 
two continua had identical shapes, with image A brighter 
than B by 0.35 mag at the time of observation. Above 
5300 A, however, the flux from B rose more steeply than 
that from A, and they were equal at ~6500 A. They also 
discovered a low-ionization absorption system at z abs = 
1.390 and noticed that in the spectrum of image B, the 
absorption lines appeared to be weaker than the corre- 
sponding lines in the spectrum of image A. This latter result 
was soon confirmed by the MMT spectra of Weymann et al. 
(1979). 

WCW suggested that differential reddening along the two 
lines of sight might be responsible for the different contin- 
uum shape in the two images. In particular, they pointed 
out that the observed break at 5300 A corresponded to a 
wavelength of 2200 A in the QSO rest frame, and 2200 A is 
a well-known resonance feature in interstellar extinction 
caused by dust grains in the Galaxy. Following this sugges- 
tion, Wills & Wills (1980) constructed various dust extinc- 
tion models to fit the observed continuum shapes. It is 
obvious that for the dust extinction model to work at 
optical wavelengths, the extinction must be stronger along 
the line of sight to image B, since the flux of image B falls 
faster at shorter optical wavelengths. 

Young et al. (1980) found a cluster of galaxies around 
0957 + 561 A, B at redshift -0.39. They attributed the 
break in the spectrum of image B to the presence of the 
brightest member galaxy, Gl, superposed on that image. 

1 Based on observations with the NASA/ESA Hubble Space T elescope , 

obtained at STScI, which is operated by AURA, Inc., under NASA con- 
tract NAS5-26555. 


Direct imaging under superb seeing conditions by Stockton 
(1980) resolved Gl. Further observations by Young et al. 
(1981a) revealed the variability of the two images. The red- 
shift of Gl was redetermined as 0.36. The evidence was 
convincing, and Wills & Wills (see their Note added in 
proof) agreed with the Gl galaxy explanation. 

We present UV spectra (2200-3300 A) of 0957 + 561 A 
and B obtained with the Faint Object Spectrograph (FOS) 
on board the Hubble Space Telescope (HST). 2 In this spec- 
tral region we found that the flux level of image A rises more 
slowly toward the short-wavlength end than does that of 
image B. We also found that the absorption system (z abs = 
1.391 1) near z em is a weak, damped Lya system with strong 
Lya absorption lines present in both images. However, it is 
clear that the two H i column densities are different, with 
the sight line to image A intersecting a large column density. 

While the difference in continuum shape in the optical 
range can be explained by contamination from the lensing 
galaxy Gl, this explanation will not work in the UV. We 
propose that the difference in the continuum shapes in the 
UV is due to differential dust grain reddening by the 
damped Lya system at z abs = 1.391 1. In this paper we fit the 
observations by assuming that dust is mixed with the gas in 
the damped Lya clouds. The direct result of this procedure 
is a determination of the dust-to-gas ratio in the damped 
Lya system. We show that by combining optical observ- 
ations with HST spectra the shape of the dust grain extinc- 
tion curve may also be determined. We compare our results 
with previous work which determined the dust-to-gas ratio 
in damped Lya clouds. 

2. OBSERVATIONS AND DATA REDUCTION 

On 1995 January 26, UV spectra (2222-3277 A) of both 
images of the gravitational lens 0957 + 561 were obtained 
with the moderate dispersion G270H grating (R ^ 1300) 
and the Red Digicon of the Hubble Space Telescope s 

2 Program 5683. 
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COSTAR — corrected Faint Object Spectrograph. The total 
exposure for image A was 6500 s and that for image B (the 
South West image with the lensing galaxy about 1" away) 
was 6560 s. Since the spectra were sampled at quarter-diode 
intervals to improve the resolution, the actual exposure 
times per pixel are one-quarter of these total exposure 
times. The observations were made through a circular aper- 
ture of diameter 0'.'43. The target aquisition sequence was 
designed to center each object to within 0'.'04 (0.5 pixel) in 
order to minimize errors due to flat-field variations. An 
additional small error in the position of the spectrum on the 
diode array may arise due to nonrepeatability in the setting 
of the filter/grating wheel (see below). 

The UV absorption lines corresponding to the two opti- 
cally identified absorption systems (see Young et al. 1981b 
and references therein) are present in both HST/FOS 
spectra. For the stronger absorption system near z em , 
Young et al. (1981b) determined an absorption redshift, 
z abs = 1.3911. We adjusted our wavelength scales to set the 
relatively weak damped Lya systems in A and B to be at 
redshift 1.3911. The shifts differ by 0.84 pixels (0.43 A). Both 
spectra are plotted in Figure 1. 

3. MODEL CALCULATIONS AND RESULTS 
3.1. Continuum Shape and Dust Grain Extinction 

Following Fall & Pei (1989) we write the optical depth 
due to dust grain extinction as 

x(X) = x B /4(A)//4(4400 A) , 

where A(A) is the extinction in magnitudes and t b is the B 
band optical depth in the rest frame of an absorber. The 


dimensionless dust-to-gas ratio, k , is defined as (Fall & Pei 
1989) 

x B = kN m / 10 21 , 

where JV HI is the H i column density in units of cm 2 . For 
the A(A) extinction curve we consider two models. Model I 
assumes a simple A(X) cct x cc 1/A law and this leads to 

/ N h | V 4400 A\ 

T(A = VKFVV — I — ) (m ° dd ' 

In model II we use the mean Galactic extinction law derived 
by Cardelli, Clayton, & Mathis (1989; hereafter CCM). 
CCM give the parameterized analytical form of A(k)/A v for 
wavelengths k > 1000 A. Since the mean diffuse ISM extinc- 
tion, A(X)jA v , above 1000 A is characterized well by the 
value R v = A V /E(B— V) = 3.1 (CCM), we adopt this value 
for Model II. We then have 


z(k) 


= k ("»i 

Vio 21 , 


\[mj Ay 

/L Ay J[_/4(4400 A)_ 


(model II) . 


In our case, k = k z = k ob J(l + z abs ), where z abs denotes the 
absorption redshift of the absorber (z abs = 1.3911 in our 
case). We have extrapolated the CCM expression below 
1000 A and applied it to shorter wavelengths near the 
Lyman limit (k = 912 A) in our calculations. The Fall & Pei 
Galactic extinction curve is very similar to our CCM curve, 
while our simple model I curve falls between the LMC and 
SMC curves in the wavelength region they used. Thus 
model I, although very simple, does approximate some real- 
istic dust grain extinction laws. 



Fig. 1. The HST/FOS spectra of Q0957 + 561 A, B. Lya lines in the strong absorption line systems are marked. The upper panel is the spectrum of 
image A, and the lower panel is the spectrum of image B. A portion of the spectrum near 2680 A which was contaminated by a noisy diode has been omitted. 
The 1 a errors are shown as dotted lines. 
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The model extinction curves are plotted in Figure 2, nor- 
malized to the B band (A B = 4400 A). The solid line is the 
CCM curve (model II) and the dashed line is model I. The 
vertical dashed line represents the short-wavelength limit of 
the CCM parameterization. 

Now let /£ bs and be the observed flux from the 
images A and B, respectively. Let /£ b ,(0) and /® o J0) be 
the intrinsic flux before absorption by dust grains in the 
damped Lya system. We have 



where N A and N B are the H i column densities along the 
lines of sight to A and B, respectively. For gravitationally 
lensed images, if there is no time delay, /£ b8 (0) and/^JO) 
have the same shape, and the first term on the right-hand 
side of the above equation is a constant, independent of 2 0bs . 
In reality, however, both time variability and time delay are 
often observed in gravitationally lensed systems. The above 
statement is still correct if the variability only changes flux 
levels and does not change spectral shapes in the wave- 
length range of interest. We assume this is the case for this 
calculation. 

In Figure 3, In [/t b ,/(0.755/^ o J] is plotted as a function 
of 4400 A/i obs . Since the original spectra of A and B are 
offset sightly in wavelength, we have interpolated to calcu- 
late flux ratios at wavelengths common to both images. The 
normalization constant, 0.755, is the flux ratio >fx oh Jf* ohs , 
estimated in the wavelength range 3000 3300 A. In this 
plot, we excluded the emission lines since they are produced 
in an extended region much larger than that of continuum 
source and may be amplified differently in images A and B. 
We also excluded the regions with strong absorption lines 
and a short region affected by an intermittent diode. We see 
clearly in Figure 3 that the flux of image A is weaker at 
shorter wavelengths relative to image B. This can be easily 



Fig. 2. — The dust grain extinction curves. The solid line is the CCM 
mean Galactic extinction curve (model II), and the dashed line is the simple 
1/2 extinction law (model I). The vertical dashed line indicates the shortest 
wavelength in the CCM determination of the mean extinction curve. The 
horizontal line shows our wavelength coverage in the absorber frame. 


^obs 
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Fig. 3. — Observed flux ratio as a function of 4400 A/2 obs . The solid 
line is the model II fit and the dotted line is the model I fit. Notice that the 
2200 A feature in the model II curve is present, but weak. In the optical 
region, model I and model II predict quite different flux ratios. 


explained as the result of differential dust grain extinction if 
the H i column density along the line of sight to A is greater 
than toward B. 

3.2. Column Density Determinations 

The values of the H i column densities, N A and JV B> can be 
measured using the damped Lya lines at z abs = 1.3911 
present in each spectrum. For each spectrum, an accurate 
estimate of the unabsorbed Lya emission-line profile on the 
blue side is important for determining the column density 
since it affects the profile determined for the damped Lya 
absorption line. We assumed that the intrinsic Lya emission 
was symmetric about z em , and reflected the right half of the 
Lya emission profile to the left half. The wavelength and 
height of the peak position of each Lya emission line were 
determined by trial and error to provide the best Voigt 
profile fit to the corresponding damped Lya absorption. 
Our results, N A = 1.9 x 10 20 cm -2 and N B = 8 x 10 19 
cm -2 , are plotted as thick solid lines in Figures 4 a and 4b. 
In Figure 4a the dashed line is for N A = 2.2 x 10 2 ° cm“ 2 
and the dotted line is for N A = 1.6 x 10 2 ° cm -2 . In Figure 
4b the dashed line is for N B = 1.0 x 10 20 cm -2 and the 
dotted line is for N B = 6 x 10 19 cm -2 . These values were 
the extremes of what would be considered reasonable fits. 
We have used an instrumental resolution FWHM = 2.0 A 
in our fits. The results are not sensitive to the resolution 
because the lines are quite wide. Our column density deter- 
minations are both higher than a measurement made in the 
low-resolution IUE spectra of the two images combined, 
]V = 7xl0 19 cm“ 2 ± 30% (Turnshek & Bohlin 1993). 

The flux level at the center of the damped Lya absorption 
lines does not go to zero in our FOS spectra. This could be 
either because the scattered light and/or dark background 
levels were underestimated or because the damped Lya 
absorbing cloud is too small to cover the whole quasar 
emission-line region. Recently Barlow et al. (1996, private 
communication) found many cases of saturated absorption 
lines with redshifts near the quasar emission redshift with 
nonzero residual flux levels at the line centers, even at the 
very high spectral resolution achieved by the Keck tele- 
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Fig. 4 .—(a) The fitted Voigt absorption line profiles for the observed (thin solid line) damped Lya absorption in the image A. The thick solid line is our best 
result N Hl = 1.9 x 10 20 cm -2 , the thick dashed line is iV HI = 2.2 x 10 20 cm -2 * and the thick dotted line is N Hl = 1.6 x 10 2 ° cm -2 . ( b ) The fitted Voigt 
absorption line profiles for the observed (thin solid line ) damped Lya absorption in the image B. The thick solid line is our best result N H j = 8 x 10 19 cm " 2 , 
the thick dashed line is N H ,= 1 x 10 20 cm' 2 , and the thick dotted line is N H j = 6 x 10 19 cm -2 . 


scope. Incomplete coverage of the quasar emission-line 
region by absorbing clouds is their favored explanation. 
Our case is similar to theirs in that z abs ~ z em , but the lines 
in 0957 + 561 have much larger column densities. In either 
case, it would be appropriate to subtract the residual fluxes 
to force the damped line centers to go to zero level. The 
residual fluxes subtracted are 3.5 x 10“ 17 ergs cm -2 s” 1 
A -1 in A and 5.8 x 10” 17 ergs cm” 2 s 1 A” 1 in B. The 
residual flux level in B sets an upper limit for the contami- 
nation by the lensing galaxy Gl. 

3.3. Results 

We have used the FIT routine in Press et al. 
(1992) to fit our model predictions to the observed flux ratio 
as a function of wavelength. Denoting a = ln{/£ b (0)/ 
[0.755/- 1 JO)]} and b = (1 - N b /N a )(/cV a / 10 2I M1 +V abs ), 
we find a = 0.181 ± 0.015 and b = 0.146 ± 0.0097 for 
model I. For model II we find a = 0.051 ± 0.0078 and 
b = 0.080 + 0.0062. These results, plus the measured iV A 
and N Bi correspond to k = 0.55 and /? . (0) =1.11 

for model I, and k = 0.31 and /? J0)//£J0)° = 1.26 for 
model II. The two best fits are plotted in Figure 3 using 
a dotted line for model I and a solid line for model II. 

We have assumed <r(y) is a constant, where y = 
In [/l b9 /(0.755/^ o J], for the FIT routine. a(y) includes the 
statistical error and the spikes produced by weak absorp- 
tion lines (we have excluded strong absorption lines) whose 
profiles in one spectrum are not exactly the same as the 
corresponding ones in the other spectrum. Judging from 
Figure 3, our assumption is reasonable. Taking estimates 
g(N a ) = 3 x 10 19 cm -2 and <t(jV b ) = 2 x 10 19 cm” 2 , our 
estimations of the error are a(k) ~ 0.18 for model I and 
a(k) ^ 0.10 for model II. Other possible sources of error 
which have not been included are uncertainty in the dust 
grain extinction curve, errors in the continuum levels used 
to measure the Lya absorption and possible weak emission 
lines in our observed wavelength region. 

In our determination of the dust-to-gas ratio, we have 
made several assumptions. We have assumed the environ- 
ments are similar and k is the same along both lines of sight 


because of the small separation at the damped Lya system. 
(At z abs = 1.3911 the separation between the two lines of 
sight corresponds to a linear size of 0.30 h$Q kpc for q 0 = 
0.5, and 0.34 h 5 Q kpc for q 0 = 0.1.) From Figure 1 we see 
that Lya absorption in the z abs = 1.1249 system (Young et 
al. 1981b) is not very strong, and we can safely neglect its 
dust grain extinction contribution. The high Galactic lati- 
tude (b ~ 48°) and small separation suggests that differen- 
tial reddening due to the Galaxy is also negligible. Some 
groups have argued that low-ionization species, such as Zn 
n and Cr ii, can give estimates of both metallicity and dust 
grain content in damped Lya clouds (see Meyer & York 
1992; Pettini et al. 1995; Wolfe 1995). With the Keck tele- 
scope it will be possible to observe these lines in both com- 
ponents and to compare k derived that way with our results. 
UV monitoring can test our assumption that the intrinsic 
continuum slopes of the two spectra are the same. 

4. COMPARISON WITH PREVIOUS WORK AND 
DISCUSSION 

Since damped Lya systems are thought to be due to the 
progenitors of galaxies like our own, they are prime candi- 
dates to contain dust grains. A direct indication that dust 
grains exist in damped Lya clouds comes from the 
reddening of background quasars (Fall, Pei, & McMahon 
1989; Pei, Fall, & Bechtold 1991 ; Fall & Pei 1995; hereafter 
Fall et al.). Fall et al. found that quasars with damped Lya 
systems in the foreground tend to appear redder, in a sta- 
tistical sense, than those without damped Lya systems in 
the foreground. Their best determinations for the dust-to- 
gas ratio in damped Lya clouds are k ~ 0.357o!o9 using the 
Galactic extinction law, /c^0.09to!o2 using the LMC 
extinction curve, and k ~ 0.06 using the SMC extinc- 
tion curve (Pei et al. 1991). In comparison, the observed 
dust-to-gas ratios in the Galaxy and the Large and Small 
Magellanic Clouds are, respectively, k = 0.79, k = 0.19, and 
k = 0.05 (see Fall & Pei 1989 and references therein). 

Before comparing the Fall et al. results with our measure- 
ments, we point out that their determinations are made in 
the wavelength region between the Lya and C iv emission 
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lines. It is interesting to see that our model II value, 
k = 0.31, is very similar to their k = 0.35 for the Galactic 
extinction law, given the fact that these values are obtained 
by two completely different methods. But our model I value 
of k = 0.55 is significantly larger than their LMC and SMC 
k values. Also our model I k is larger than in model II while 
Fall et al. get smaller k values for LMC/SMC type extinc- 
tion than Galactic type extinction. The main reason for the 
disagreement is that we have a smaller intrinsic flux ratio, 
/lb,(0)//£>b,(0), in model I (1.11) than in model II (1.26), and 
to produce the observed f\Jf\ ob> ^ 1.32, we need a larger k 
for a model I cloud. 

Now let us examine why our model I result is so different 
from the Fall et al. LMC/SMC results. One possibility is 
that the LMC/SMC type extinction curve is quite different 
from the simple 1/A law at wavelengths shortward of Lya 
emission. If this is true, then our model I result will have to 
be recalculated. At present we know little about the 
LMC/SMC type extinction curve in the FUV and EUV 
regions and look forward to future observations that may 
improve the situation. 

Another possible explanation for the discrepancy 
between our large model I k value and the small k value of 
Fall et al. is that there is a large dispersion in k. Fall et al. 
assumed that there is a single k for all damped Lya 
absorbers because it is difficult to model the dispersion in k. 
However, a large dispersion in k is expected for at least two 
reasons: (i) damped Lya absorbers are found in various 
states of chemical evolution; (ii) radial abundance gradients 
are probably present, and the scatter in impact parameters 
implies a large spread of k. Boisse (1995) argued that the low 
dust-to-gas ratio in damped Lya absorbers may be the 
result of a selection effect; a smaller mean k value would be 
favored because QSOs obscured by damped absorbers with 
large k values may be too faint to be included in spectro- 
scopic samples. 

A large dispersion in k seems to be consistent with the 
observations of metal abundances in damped Lya systems 
by Pettini et al. (1995 and references therein). They used Zn 
as a tracer of metallicity since Zn is relatively undepleted by 
dust grains, and Zn should also trace Fe, the most widely 
used indicator of metallicity (Pettini et al. 1995; Wolfe 
1995). Pettini et al. found that, while on average the metal 
abundance in the damped Lya systems at redshifts z ~ 2-3 
is about 1/10 of the solar value, there is a large dispersion, 
from near solar to less than 10" 2 of the solar value, in the 
degree of metal enrichment attained by different damped 
absorbers at essentially the same epoch. These results are 
confirmed by the recent observations of Wolfe & Prochaska 
(1995, private communication) using the HIRES echelle 
spectrograph on the Keck 10 m telescope. Although it is not 
known how the metal abundance relates to k directly, it 
may be that a large dispersion in metallicity means a large 
dispersion in k. So at the present time we cannot rule out 
the large model I k value for the 0957 + 561 A, B absorption 
system. Measurements of Zn and Cr in other damped Lya 
systems at redshifts similar to (or lower) than that of 
0957 + 561 imply very low dust-to-gas ratios (Meyer & 


York 1992; Steidel et al. 1995; Meyer, Lanzetta, & Wolfe 
1995). If the value of k we measure here is correct, this 
implies that there is a considerable range of k at low red- 
shifts as well. A measurement of Zn and Cr in 0957 + 561 
would be an important test. 

We should also point out that this damped system red- 
shift is near the quasar emission redshift. It has been sug- 
gested that metal absorption systems at z abs ~ z em could be 
associated with the QSO itself either because they arise in 
galaxies hosting the QSO or because they are the signature 
of ejected material (see Petitjean 1995 and the references 
therein). By studying nine absorption systems in two QSOs, 
0424—131 and 0450—131, Petitjean & Bergeron (1994) 
found that the systems within 10,000 km s _1 of the QSO 
have abundances similar to or even in excess of the solar 
value. In 0957 + 561A, B, the difference between z abs and z em 
corresponds to about 3000 km s" and the large model I k 
value might indicate an enhanced metal abundance near the 
quasar. 

Fall et al. and many others (see the references cited in the 
Fall et al. papers) have tried but failed to detect the 2200 A 
feature of Galactic type dust in damped Lya absorbers. In 
Figure 3 we extend the fits for both models into the optical 
to show this region. It is clear that for k = 0.31, observ- 
ations with exceptionally good S/N would be required to 
detect the 2200 A feature of Galactic type dust grains. 
However, the two models predict quite different flux ratios 
in the optical region, and it should be possible to select 
between these models for 0957 + 561. Observations with the 
HST can exclude any contribution from Gl, which cur- 
rently contaminates the optical spectrum of 0957 + 561 B, 
while obtaining the entire optical/UV spectrum simulta- 
neously. If, in addition, observations of A and B are separat- 
ed by the time delay (see, e.g., Kundic et al. 1995), it will be 
possible to eliminate errors due to any differences in contin- 
uum slopes caused by spectral variability. 

Since it may strongly affect our view of the distant uni- 
verse, dust grain extinction has important cosmological 
consequences. This topic has been studied by several 
groups. In these studies dust grains are either associated 
with intervening galaxies (Ostriker & Heisler 1984; Heisler 
& Ostriker 1988; Wright 1986, 1990) or attached to quasar 
absorption-line systems (Fall & Pei 1995; Zuo & Phinney 
1993; Ostriker & Vogeley 1990). The advantage of using 
quasar absorption-line systems is that these are observed in 
large quantities at high redshifts. The outcome of the 
obscuration calculations depends critically on the k value 
used. When more damped Lya absorbers in gravitationally 
lensed systems are discovered, the method presented here 
will become a powerful tool for determining the value and 
dispersion of /c, as well as the shapes of extinction curves, 
especially in the FUV and EUV regions. 
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